Parameterized Post-Friedmann Signatures of Acceleration in the CMB 
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We extend the covariant, parametrized post-Friedmann (PPF) treatment of cosmic acceleration 
from modified gravity to an arbitrary admixture of matter, radiation, relativistic components and 
spatial curvature. This generalization facilitates the adaptation of Einstein-Boltzmann codes for 
solving CMB and matter perturbations in the linear regime. We use such a code to study the 
effect of metric evolution on the CMB through the integrated Sachs- Wolfe effect. We discuss the 
ability of modified gravity to alter the low multipole spectrum, including lowering the power in the 
quadrupole. From a principal component description of the primary metric ratio parameter, we 
obtain general constraints from WMAP on modified gravity models of the acceleration. 

PACS numbers: 



I. INTRODUCTION 

In the absence of theoretically compelling models for 
the acceleration of the expansion, it is useful to have a 
parametrized description of possible deviations from the 
standard cosmological constant model. To explore mod- 
ified gravity explanations of the acceleration we need a 
parametrized post-Friedmann (PPF) description of grav- 
ity that parallels the parametrized post-Newtonian de- 
scription of solar system tests. 

Many attempts toparametrize deviations phenomeno- 
logically (e.g. [3,3, 4,[a|) or across a limited range of 
scales (e.g. [Q, 0, exist in the literature. Hu 

& Sawicki [ll| recently introduced a PPF approach that 
spans all linear scales and includes an ansatz for non- 
linear phenomenology. The linear framework is based 
on enforcing the requirements of a metric theory, with 
small deviations from the Friedmann-Robertson- Walker 
metric, and covariant conservation of the stress energy 
tensor. This approach provides an excellent description 
of at least two specific models of modified gravity, the 
self- acceleratiiig branch of the braneworld acceleration 
model (PGP |13 |) and the modified action f{R) model 

In this Paper, we extend the PPF approach to in- 
clude relativistic matter components and spatial curva- 
ture. This generalization yields a self-consistent means of 
evolving metric perturbations from initial conditions in 
the radiation dominated era through to the present. We 
also phrase the PPF description in a manner that is useful 
for adapting an Einstein-Boltzmann code (e.g. [H, [T7| ) 
for modified gravity. Furthermore, specific modified 
gravity models already require such components at low 
redshift, e.g. the self-accelerating branch of DGP re- 
quires spatial curvature to fit distance measures [l^ and 
the normal branch requires dark energy to explain accel- 
eration. 



We use this generalization to study the impact of evo- 
lution in the metric on the cosmic microwave background 
(CMB) through the integrated Sachs- Wolfe (ISW) effect. 
In the DGP and f{R) models of cosmic acceleration it is 
well known that the ISW effect provides one of the most 
powerful and robust tests available @, [H [H, HO, [SI, HI . 
Since it arises from fluctuations near the horizon scale 
during the acceleration epoch, it requires a fully covari- 
ant description such as the one developed here. It also 
therefore tests gravity on the largest scales observable. 

We begin in fjllj with the multicomponent generaliza- 
tion of the PPF description. This derivation is given in 
the comoving gauge but we describe in the Appendix the 
corresponding relations in the synchronous gauge. In mill 
we explore the phenomenology of the ISW effect in the 
CMB and examine current constraints. We discuss these 
results in TOdI 



II. MULTICOMPONENT FORMALISM 

In mi A[ we adapt the structure of multicomponent 
covariant linear perturbation theory to modified gravity 
models. This treatment reduces to that of 11] in the 
limit of a flat, matter only universe. We discuss gauges 
relevant for the PPF parametrization in ^11 Bl and derive 
the explicit representation in mi CI 



A. Covariant Structure 

For the purposes of constructing a self-consistent co- 
variant parametrization of modified gravity, it is conve- 
nient to view the modifications to gravity in terms of an 
additional "dark energy" stress tensor 0, [2^ . Given a 
metric theory of gravity, we are always free to define the 
dark energy stress tensor to be 
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where we use the subscript "T" for the total stress energy 
tensor, which combines all true components, and "e" for 
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the effective dark energy here and below. Note that we 
allow the total stress energy tensor to contain true dark 
energy components such as a cosmological constant or 
scalar fields for maximal generality (c/. 

The effective dark energy behaves as a separately con- 
served system by virtue of the Bianchi identities and the 



conservation of T!^'^ [24 



^ " 8^G 



V^T^" = . 



(2) 



However the closure relation for this conservation law, 
which specifies the relationship between the components 
of the stress energy tensor, can and will in general depend 
on the matter content in contrast to minimally coupled 
dark energy models ful,!!!!!!- 

As usual, statistical isotropy requires that at the back- 
ground level the stress tensors can be parametrized by 
the energy density and pressure 



(3) 



for both the total and the effective stress energy tensor. 

The background Einstein tensor is built out of the 
Friedmann- Robertson- Walker metric 

ds"^ = gfj.vdx'^dx" = a^{-drf' + ^^jcix^'dx^) , (4) 

where 77 = J dt/ a is the conformal time. In spherical co- 
ordinates 7ij , the spatial metric with constant curvature 
K, can be represented as 



-/.jdx'dx^ = dD^ + D\dn , 



(5) 



where Da = K~-^/^ sm{K^/^D) is the angular diameter 
distance. The Einstein equation ([1]) becomes the usual 
Friedmann equation 



The conservation laws 



Pt 



-3{pT +Pt) ■ 

-3{pe+Pe) , 



(6) 



(7) 



close the system of equations for the background equa- 
tions. The effective dark energy is thus parametrized in 
the same way as true dark energy: by the density today 
in units of critical fie = 87rGpe(lna = 0)/3Hq and the 
effective equation of state Weihia) = Pe/Pe- 

Scalar linear perturbations may be decomposed into 
the eigenfunctions of the Laplace operator 



and its covariant derivatives 

= (-fc)v,r, 



Y 



(fc- 



(8) 



(9) 



The most general scalar linear perturbations to the met- 
ric of a wavenumber k can be parametrized by j26l . \2l\ 

Sgoo = -a^{2AY), 

Sgot = -a^BYi , 

Sg,, = {2HlYj,, + 2HtY,, ). (10) 

Likewise the stress energy tensors can be parametrized 
as 

5T\ = -5p, 

ST,' ^~{p + p)vY\ 

dT)^6pYd'j+pIlY). (11) 

We again describe the modification to gravity with an ef- 
fective dark energy stress tensor and allow the total stress 
energy tensor to be composed of multiple components 

6pT 



[pi 



Pt)vt 
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(12) 



By definition, Eqn. ([T|) enforces the usual 4 Einstein field 
equations 28] 



3 kn kjj 



A + Hl 
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H^CKkl 

+ - + 

BttG 



6p + i{p + p)- 



B 



B' + 2B 








kn 






'^H J 



H-^kjj 

H'rp 



K 



AttG 



{p+p) 



[aHf 
V- B 



B 

kH 



A' 



2 + 2 i^] A - —(B' + B) 

H 3 3 ^ ' 



1 



H';,-i2 + ^)H', = ^iSp+-5p) 



(13) 



where the stress energy components on the rhs combine 
the total and effective contributions. Here ' = d/d\na, 
kn = {k/aH), and ck ~ ^ — iK/k^. The conservation 
laws become the continuity and Navier-Stokes equations 



dp'i + 3{5pi + 5p,) = ~{p, + Pi){kHVi + 3H'^) , 



(14) 



[a\p, + p,)iv, ~ B)y 
a^kn 



= Spi - -CKPi^i + {pi +Pi)A , 



for all components separately conserved including the ef- 
fective dark energy. 
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B. Gauge 

The covariant Einstein and conservation equations ([Tc 
and (fH|) apply to any choice of gauge. For the PPF con- 
struction it is useful to work with variables that take on 
certain meanings such as curvature and potential fluctu- 
ations in specific gauges. For numerical codes it is useful 
to have a covariant representation of such variables so 
that they may be accessed from other gauge choices. 

Under a gauge transformation defined by the change in 
conformal time slicing T and spatial threading L ^26j, |27i] 

+ (15) 
= + LY' , 

the metric variables transform as 

A = A- aH{T' + T) , 
B = B + aH{L' + knT) , 

Ht = Ht + aHkuL , (16) 
and the stress energy components transform as 

5p ~ 5p — p'aHT , 

dp ~ 5p — p'aHT , 
V = V + aHL' , 

n = n. (17) 

A gauge is fully specified if the functions T and L are 
uniquely defined. 

Following fll'l, we shall construct the PPF descrip- 
tion from a combination of (total matter or) comoving 
and Newtonian gauge quantities. The comoving gauge is 
specified by the conditions 



B = VT , 

Ht = 0. 



(18) 



They fully specify the gauge transformation from an al- 
ternate gauge choice 



T = (wT - B)/k, 
L = -Hr/k . 



(19) 



To avoid confusion between fluctuations defined in dif- 
ferent gauges, we will assign special variable names to 
comoving gauge quantities 



c = 


Hl 




A. 




Sp, 


Ap = 


Sp, 


V = 


V . 



(20) 



Ap should not be confused with pA = p{Sp/p). 

Similarly, the Newtonian gauge is defined by the con- 
dition B = Ht = and the transformation 



T 



L 



B 



H'rp 



kk 



H 



(21) 



To avoid confusion we define 



Hl- 
A. 



(22) 



The relationships between the two metric fluctuations are 

(23) 
(24) 



kn 



k 



Vt 



H 



We refrain from utilizing matter variables in Newtonian 
gauge but note that velocities in the two gauges are the 
same. We discuss the synchronous gauge representation 
in the Appendix. 



C. PPF Parameterization 

With the full covariant framework of linear perturba- 
tion theory in place, we now generalize the PPF descrip- 



tion of the effective dark energy stress tensor [11[ for 
multiple relativistic components and spatial curvature. 

Following the original construction, we demand that 
the additional PPF contribution satisfy two require- 
ments. On superhorizon scales kn = k/aH <ti 1, the 
curvature ^ in the comoving gauge is only altered by the 
effective dark energy at second order in kn- Conserva- 
tion of C in the absence of curvature and non-adiabatic 
stress fluctuations is a consequence of energy-momentum 
conservation 12811 and applies to modified gravity models 
that satisfy it [2g |. 

The third Einstein equation reads 



K Vt 
{aH? k^ 



AttG 



ifie+Pe) 



V,-Vt 



(25) 



and the Navier-Stokes equation for the total velocity ([T 
gives 



i = - 



Apj 



ckPt^t 



Pt +Pt 



(26) 



Without further loss of generality we can parametrize 
the effective dark energy contribution in this limit by a 
function f(^{a) where 

AttG , Ve-VT 1 , . ^, 
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since Vr = 0{kHC) for adiabatic fluctuations. The equa- 
tion of motion for C. in this hmit now reads 



lim Q = —kuVT 



fcjf<i 



PT +PT 



fc2 



(28) 



The second condition is that the metric satisfies a 
Poisson-like equation in the kn 3> 1 quasistatic limit 

47rG AtPt + cxPt^t 



lim 

kH-»i CKkjjH^ l + fcia) 



(29) 



where fa depends on time alone and $_ = ($ — 'I')/2. 
To make these two limits compatible wc introduce a pa- 
rameter r such that 



47rG 



CKk%H^ 



[AtPt + ckPt^t] (30) 



on all (linear) scales. Comparison with the first Ein- 
stein equation for the Newtonian metric perturba- 
tions gives the first closure relation for the effective dark 
energy 



PeAe + 3{pe+Pe) 



Ve-Vr 



- CKPe^e 



k'^CK 

47rGa2 



r. (31) 



Given that modified gravity models are most simply 
parametrized in terms of the relationship between $ and 
that they induce, we describe the second closure con- 
dition through the effective anisotropic stress 



(32) 



H 



Parametrizing this relation with a free function g{a, k) 
again loses no further generality (c/. 0]). Without a 
modification to gravity g = Q. If the true anisotropic 
stress is negligible then g = and for this reason 

we will refer to it as the metric ratio parameter. Equa- 
tion (15^ defines the effective anisotropic stress as 



4^G 

rrPelle 



mk 



(33) 



H 



To complete these equations we must determine an 
equation of motion for F that is consistent with the 
two requirements ([28]) and ([29|) . Taking the derivatives 
of Eqn. (1221) and ([50]) with the help of the continuity 
and Navier-Stokes equations, we obtain the condition at 
A;// < 1 



lim r' = s* - r . 



(34) 



where the source 



47rG 



r{ff[(PT7rT)' 



5 + 1 {g + l)klH^ 

' [{9 + /c + 9fc){PT +Pt) — {Pe + Pe)] kn Vt] 



(35) 



Note that S* = if the modification to gravity vanishes: 
5 = and = 0. 

In the opposite limit a comparison of Eqns. ([29]) and 
imply 



lim r = /G$_ 



(36) 



To satisfy both limits, we take the equation of motion for 
r to be 

{l + 4kl)[V' + T + clkl{T- fG-^^)] = S. (37) 

Substituting this relation into the derivative of the Pois- 
son equation, we obtain 



where 



ApT — ^ckPt^t K 



PT +PT 



3 + 1 
F 



5 - F' - F + /c 



kHVT (38) 
ATrG{pT + Pt) Vt 



F(a) 



1 



-^(5 + l)-?:^ [PT+PT) . 

K Ck 



(39) 



Comparison with Eqn. (|25p implicitly defines the effective 
dark energy momentum density 



Ve^Vr 
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1 



k 



H 



47rG(pe - 

S-T' 



■Pe) F 



(40) 



47rG(/9T + Pt) Vt 



Likewise the effective pressure can be defined from this 
relation through the dark energy Navier-Stokes equation 



[a\pe+Pe)[Ve-VT)]' 

a^kn 



Ape - -CKPe^e + {Pe + Pe)C ■ 

(41) 

This completes the multicomponent generalization of the 
PPF parametrization. It is described by three free func- 
tions /(^(Ina), fcllaa) and g(liia,k) and one parameter 
cr in addition to the usual We{a) — Pe{a) / pe{a) which de- 
termines the background expansion. In the limit that the 
universe is spatially flat and relativistic matter compo- 
nents are subdominant, this description exactly matches 
the original treatment of [ll| . 

In our numerical implementation, we evolve the density 
Ai and relative velocity Vi — Vt components through the 
comoving gauge conservation equations and the Boltz- 
mann hierarchy for the radiation closure condition and 
anisotropy using the code of [13, [3l|. Specifically we 
evolve the comoving curvature fluctuation through 
Eqn. ([55)1 and ^ from the constraint Eqn. (j^ . The final 
metric perturbation is the total velocity Vt(= B) which 
is evolved using the combined Newtonian and comoving 
Navier-Stokes equations [see also the gauge relation 



V^ + VT = kH{-^-0: 



(42) 
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FIG. 1: CMB temperature power spectrum as a function of 
the amplitude of the metric ratio parameter today go given 
the evolution described by Eqn. (|44|) . ACDM corresponds to 
go — 0. Increasing \go \ in the positive direction monotonically 
increases the ISW effect at low multipoles whereas in the neg- 
ative direction it first decreases and then increases the effect. 



The other parameters have been fixed at 
and We = —1. 



0.01, 



cr 



where $ and ^ are specified by Eqn. ((30| and g{\na, k). 
Finally F is obtained by evolving Eqn. ([57)1 . The second 
gauge relation of Eqn. between ( and $ is used as a 
test of numerical accuracy. 

Note that the relationship between F and g{lna,k) 
and the effective energy density, anisotropic stress, mo- 
mentum density and pressure fluctuations given by 
Eqns. dni), (ESI), (SOI) and dU) are only imphcitly used. 
Likewise while the effective stress energy obeys the con- 
tinuity and Navier-Stokes equations, they are not used in 
the numerical scheme. These relations are however useful 
for adapting the numerical scheme for the synchronous 
and other gauges (see Appendix). 



III. PPF ISW EFFECT IN THE CMB 

In general, our PPF description of cosmological de- 
viations from general relativity includes one function of 
space and time that interrelate the metric fluctuations 
g{\na,k), two functions of time only fc (In a) and (In a) 
that determine the relationship between the metric and 
the matter, and one parametrized scale cr that bridges 
the transition between the latter two relations. These 
free parameters allow a complex range of phenomena 
for the scale-dependent evolution of linear perturbations. 
The impact of such deviations on specific observables 
however limits the relevant parameters to a more man- 
ageable range. 

Here we shall focus on the CMB temperature 
anisotropy induced by the evolution of the metric dur- 



ing the acceleration epoch, the so-called the integrated 
Sachs- Wolfe (ISW) effect. This effect also highlights the 
extensions to the PPF description introduced in the pre- 
vious sections. They allow an Einstein-Boltzmann code 
to self-consistently calculate the net CMB anisotropy 
from the initial epoch where radiation dominates through 
to the present. 

Though the formalism allows for high redshift modifi- 
cations of gravity, we work in the context that they only 
appear during the recent acceleration epoch. Hence the 
well-tested high multipole structure and polarization of 
the CMB are left unchanged. 

In ijlll Al we describe the impact of the PPF parameters 
on the low multipole temperature anisotropy of the ISW 
effect. We examine constraints on these parameters from 
WMAP in mil Bl and the reverse-engineering of models 
to produce specific features in the spectrum in mil CI 



A. Phenomenology 

We begin by specializing the PPF description for near 
horizon scale perturbations. The ISW effect is associated 
with the evolution of the metric potential <&_ during the 
acceleration epoch of z < 1 and wavenumbers k ~ 10"'^ 
Mpc"^. Correspondingly, it is important to characterize 
the PPF metric ratio parameter g{hia, k) at such epochs 
and scales. Since g{\n a, k) should become independent 
of scale once kn ^ I let us take the functional form 



ff(lna, k) 



gsH(lna) 
1 + {cgkn) 



(43) 



We shall use the parameter Cg to explore the range of 
scales that the ISW effect tests. We take 5 ^ for 
kn '3> 1 and /g = so as to restore general relativity 
on small scales. Such models would evade all current 
constraints from large-scale structure and are uniquely 
probed by the CMB. On the other hand, we shall see 
that if Cg < 0.01 our constraints are valid for other choices 
since the high k behavior of g does not impact the ISW 
effect. In general, our treatment should be interpreted as 
constraining the average g during the acceleration epoch 
and for scales near the horizon regardless of its behavior 
on small scales. 

Next, we assume that deviations in g will only appear 
when the modification to gravity becomes an important 
contributor to the expansion rate H cx p^/^. We therefore 
take a baseline functional form of 



5SH(lna) = go[ — — 



1/2 



(44) 



such that gsH(O) — go- In the self-accelerating branch of 
DGP, g grows to go 1 by the present on scales near 
the horizon and in the f{R) models g = —1/3 on scales 
below a Compton wavelength that increases with time 
[Tl| . In the next section we will explore variations from 
the baseline functional form of Eqn. 
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FIG. 2: Effect of changing the scale at which metric ratio 
deviations occur through Cg. Once Cg > 0.1 deviations have 
been suppressed for k > 10~^Mpc~^ where the ISW effect 
at the lowest multipoles peak. Other parameters have been 
fixed to go = 1, cr = 1 and We = —1. 



FIG. 3: Effect of changing the scale below which the qua- 
sistatic Poisson equation (|29p holds through cr. At a fixed 
metric ratio g, the ISW effect is enhanced by lowering cr and 
delaying the onset of the quasistatic dynamics. The other pa- 
rameters have been fixed to go — 1, Cg = 0.01 and We = —1. 



The remaining important parameter is cr vifhich de- 
fines in units of the Hubble scale where the transition 
to the quasi-static Poisson equation occurs. We begin 
by choosing cr = 1 as motivated by the DGP and f{R) 
models "IT]. 

Finally for definiteness we take /^(Ina) — OAgsii{\na) 
for the relationship between the density A and metric 
fluctuations $_ on superhorizon scales (see [ill] for a 
discussion). The ISW effect is insensitive to variations 
in this choice if cr ~ 1. 

We take We = —1 to illustrate effects that are com- 
ing purely from the modification of gravity and not the 
change in the expansion history. For the other cosmo- 
logical parameters, we take the maximum likelihood val- 
ues for WMAP3 f^n^/i^ = 0.128, nbh^ = 0.0223, 
He = 0.76, T = 0.092 with an initial power spectrum of 

f^..;( " . )"" (45, 

27r2 « V0.05Mpc"V 

with = 4.56 X 10-^ and n = 0.958. 

We begin by exploring the impact of the amplitude 
of the metric deviation go in Fig. [1] with the other pa- 
rameters set to their fiducial values: Cg = 0.01, cr = 1, 
We = —1. Note that go = is exactly equivalent to the 
ACDM model. For 50 > (solid Hues), the ISW contribu- 
tions increase monotonically since raising go corresponds 
to increasing the decay of the potential For 50 < 0, 
the decay is slowed and eventually turns into growth. 
Correspondingly, the ISW effect is first reduced and then 
enhanced as \go\ increases. 

In Fig. [21 we illustrate the effect of suppresing g at 
high k through Cg in Eqn. for gQ = 1, cr = 1 and 
We = — 1. As Cg is raised beyond 0.1, the enhancement to 



the ISW effect begins to go away. This confirms that the 
relevant scales at which g is constrained from the effect 
is near a tenth of the horizon scale or k ^ 10~^ Mpc~^. 

Changing the quasistatic transition scale through cr 
has a substantial impact on the relationship between the 
metric ratio g and the ISW effect. As discussed in jTl| . 
the superhorizon impact of g implied by the conserva- 
tion of the comoving curvature ^ tends to be substantially 
larger than the quasistatic limit would imply. In Fig. [3] we 
show the effect of decreasing cr for go = 1, Cg = 0.01 and 
We — —1. This extends the region where the superhori- 
zon behavior hold. Once cr ^ 0.1, wavenumbers relevant 
to the ISW effect are impacted leading to a substantial 
enhancement of the effect for the same g. The amount 
of enhancement then depends also on the superhorizon 
parameter and we take = OAggn throughout as 
an illustrative example. For the ISW effect f( is largely 
degenerate with cr since they both control the interpola- 
tion between the superhorizon behavior and quasistatic 
behavior both of which only depend on g [Tl| . 

Finally, in Fig.|3]we show the effect of varying the back- 
ground expansion history for a fixed go = 1, Cg = 0.01 
and Cr — 1- Here fie is adjusted to keep the distance 
to recombination fixed for the same flmh'^ and ftbh"^ of 
the fiducial model. For the relatively small variations 
allowed by cosmological distance tests, the changes are 
small compared with cosmic variance and detectable vari- 
ations in g. 



B. Constraints 

In the previous section, we have shown that for a given 
transition scale cr, the ISW effect constrains the metric 
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FIG. 4: EflFect of changing the expansion history through 
is small compared with cosmic variance and the effect of the 
PPF parameters for small deviations —1.1 < We < —0.9. The 
other parameters have been fixed to go = 1, Cg = 0.01 and 
cr = 1. 



FIG. 5: WMAP likelihood of models as a function of go for 
various cr with Cg = 0.01 and other parameters fixed. Taking 
Cr = 1 provides conservative constraints. Normalization is 
arbitrary. 



ratio parameter g ai k ^ 10^^ Mpc^^ during the ac- 
celeration epoch. We can now assess how well WMAP 
constrains the amplitude of deviations in t;. In lieu of 
a joint analysis of cosmological data sets that determine 
the distance redshift relation, we have here fixed the ex- 
pansion history and other cosmological parameters to 
ACDM. Given how well most parameters are currently 
fixed, the most important caveat introduced by this as- 
sumption is that the initial power spectrum is taken to 
be of the power law form given by Eqn. (|^5)) . We shall 
return to this point below. 

We begin by examining constraints on the amplitude 
parameter 50 .9(lna = 0) of Eqn. (I44p . The relative 
WMAP likelihood is shown in Fig.[5]for several choices of 
Cr and Cg < 0.01. Constraints tighten as cr decreases due 
to the enhancement in the ISW effect shown in Fig. [31 A 
conservative interpretation of the constraint that is also 
well-motivated by the DGP and f{R) examples would 
be to take cr = 1. The distributions are all consistent 
with (70 = 0. The statistically insignificant preference for 
negative go is associated with a slight lowering of the low 
£ multipoles favored by the data as we discuss in the next 
section. 

Even taking the conservative cr = 1 limit, there is a 
problem in interpreting the constraint in a general con- 
text. The parameter represents the metric ratio today 
and is not directly constrained by the data. We have fur- 
thermore assumed an ad hoc evolution of g in Eqn. (I44p 
motivated only by a rough scaling with the impact on 
the expansion rate. 

In order to place more general constraints on the met- 
ric ratio 5, we use a principal components construction to 
determine the temporal weights that the ISW effect ac- 
tually constrain. The principal component construction 



begins with the calculation of the Fisher matrix 

d\nCid\nCi 



dpi dpj 



(46) 



where Pi are a set of Np parameters that represent the 
values of a spline interpolated function p(lna) 



gSH(lna) =p{\na) [ 



(47) 



evaluated at equally spaced intervals in In a. The deriva- 
tives are evaluated at the ACDM model of pi = 0. We 
further take Cg < 0.01 and cr = 1 here. 

We then decompose the Fisher matrix into principal 
components indexed by fi 



(48) 



where a?, is the Fisher estimate of the variance of mode 
associated with the linear combination ^iSi^pi. Rank 
ordered in increasing variance, we find that the second 
mode has 5 times the variance of the first mode and so 
we shall use only the first mode in the following analysis. 

More specifically, let us define a continuous function 
Vi which is the high Np limit of the discrete eigenmode 



Vi{\Yiai) oc Sii , 



normalized so that 



dlnaVf 



1. 



(49) 



(50) 



The metric ratio g that this mode represents is given by 



ffi(lna)=(^^)'^Vi(lna) 

\PT Vie J 



(51) 
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FIG. 6: First principal component gi of the metric ratio pa- 
rameter as a function of scale factor a. Constraints derived 
from this component represent a temporal average of g of this 
functional form through Eqn. (|54p. Note that the maximum 
value gmax — 0.5 aX z = 0.68. 



FIG. 7: WMAP constraints on ges from the first principal 
component. Note that g^s is the maximum value of g attained 
by the principal component. For comparison the direct con- 
straint on go with cr = 0.1 from Fig. [5] is overplotted with 
the conversion geff = 0.54go from Eqn. (|54p . 



and is shown in Fig. [6l Note that the maximum value 
attained in this function is ^max = 0.5 at a redshift of 
^max = 0.68. The ISW effect thus constrains a weighted 
average of g around Zmax- The falloff in sensitivity near 
the present epoch is due to projection effects. 

We now use this functional form of g to evaluate the 
WMAP constraint on the amplitude g^s 



gsH(lna) = .gcff 



gi(lnQ) 



(52) 



Fig. [7] shows that the posterior probability distribution 
of gcff is nearly Gaussian and that 



ffoff 



-0.12 ±0.27. 



(53) 



Note that gcff is the value of gsH at the Zmax- 

Given a specific model for g, this constraint should 
be interpreted as one on a given weighted average of g 
around Zmax; specifically 



5cff = 5max / (ilnap(lna)Vi(lna) 



(54) 



For example, employing the model of Eqn. (|44p we ob- 
tain (Jeff — 0.54go- In Fig- El we show that the direct 
constraints on go compare well with those inferred from 

5eff- 



C. Designer Models and Caveats 

The constraints obtained in the previous section are 
fairly general in that a model that violates them while 
maintaining the underlying assumptions of a smooth ini- 
tial power spectrum and an expansion history close to 



ACDM will be disfavored by the data. However there 
remains a wide range of interesting phenomenology as- 
sociated with acceptable models and a weakening of the 
underlying assumptions. 

As an example of phenomenology that PPF models 
can produce that is not possible in ACDM, consider the 
statistical curiosity of the low quadrupolc power in the 
CMB. Figure [8] shows that in the maximum likelihood 
ACDM model, the observed quadrupole is near the 95% 
cosmic variance limit. A quadrupole realization as ex- 
treme as the WMAP data only occurs in ~ 5.5% of the 
time. Of course in the first say 100 multipoles one ex- 
pects and finds a few such events. The quadrupole is 
only special because it represents the largest observable 
scale and may hint at new physics beyond the standard 
cosmological model. 

Construction of a model that sharply lowers the power 
in the quadrupole but not the octopole is hindered by 
the fact that in ACDM it gets nearly equal contributions 
from the ISW effect and the Sachs- Wolfe effect. These 
contributions are separated by a decade in physical scale 
making it difficult to invoke initial conditions as an ex- 
planation of the low quadrupole [s^. Furthermore for 
dynamical dark energy models based on scalar fields, it 
is difficult to reduce the ISW effect while maintaining an 
expansion history close to ACDM [sj, [s^ . 

As we have seen, a modified gravity model can alter 
the evolution of the metric even for an expansion history 
that is indistinguishable from ACDM. Let us then con- 
struct a model which nearly eliminates the ISW effect at 
low multipoles while leaving the spectrum at > 10 un- 
changed. To enhance the impact on the low multipoles 
let us take cr = 0.1, Cg < 0.01 and take a very steep 
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FIG. 8: Temperature power spectrum of the ACDM model 
with the WMAP data and noise error bars overplotted. 
Shaded regions represent the 68% and 95% cosmic variance 
intervals for the model. A PPF model constructed to elimi- 
nate the ISW effect can lower the power in the quadrupole. 
Shown is the model discussed in the text with and without 
an initial power spectrum cut off at fcmin ~ 2.4 x fO~*Mpc~^. 



evolution of g such that deviations appear only recently 



5SH(lna) = A 



Pe 
PT 



(55) 



We find that A = —0.7 causes a growth at late times 
that nearly offsets the earlier decay leaving the ISW ef- 
fect nearly absent in the quadrupole. Because of the 
steep evolution, this model evades the constraints from 
the previous section. We show the spectrum of this model 
in Fig. [5] (dashed line). This model is a better fit to the 
WMAP data with a A21nL = -1.5. Of course given the 
number of parameters and ad hoc assumptions of this 
model, such an improvement is not statistically signifi- 
cant. On the other hand, explorations of this sort can 
suggest what types of modified gravity models would be 
fruitful to study: in this case, models with g < on the 
horizon scale near the present. 

With the ISW effect nearly eliminated the quadrupole 
can be further lowered by a change in the initial power 
spectrum of Eqn. ([45|) . Placing a cutoff of fcmin = 2.4 x 
lO"'^ Mpc^^ produces the spectrum shown in Fig. [5] and 
an improved likelihood of A21nL = —2.8. 

This example also highlights the primary caveat of 
applying the constraint on g^s in the previous section. 
That constraint assumes a featureless initial power spec- 
trum where the constraints at high i nearly fix the initial 
power at the horizon scale today. While this is reason- 
able given the expectation that gravity is only modified 
at late times, it remains an assumption nonetheless. 



D. Discussion 

We have generalized the parametrized post-Friedmann 
(PPF) description of cosmic acceleration from modified 
gravity to include multiple relativistic matter species 
and spatial curvature. This generalization facilitates 
the adaptation of Einstein-Boltzmann codes for modi- 
fied gravity. We have adapted a comoving gauge code to 
study constraints on deviations from general relativity on 
the largest scales through the CMB. 

The integrated Sachs-Wolfe (ISW) effect probes the 
evolution of the metric <&_ on scales near the horizon dur- 
ing the acceleration epoch. Under the PPF parametriza- 
tion, modifications to gravity can change the amplitude 
and shape of the ISW contributions to the low multipole 
CMB temperature power spectrum even if the expan- 
sion history is indistinguishable from that of a cosmo- 
logical constant. For example in some regions of param- 
eter space, the ISW effect can be nearly eliminated at 
the quadrupole bringing the predicted ensemble average 
quadrupole nearer to its observed value on our sky. 

We use the WMAP data to constrain such modifica- 
tions and find that deviations in the metric ratio param- 
eter g are constrained at the level of gcS = —0.12 ± 0.27 
at z ~ 0.7 and k ~ 10"'^ Mpc~-^. Correspondingly, spe- 
cific models such as the DGP and some f{R) models are 
constrained by the data at levels consistent with the find- 
ings of previous work (e.g. [H, H^l). By phrasing these 
and other constraints in the model-independent PPF de- 
scription, we gain a more general understanding of what 
aspects of general relativity are tested by the observa- 
tions. 

The modified gravity framework studied here should 
also enable future studies of other cosmological observ- 
ables on large scales that are beyond the scope of this 
work such as the lensing of the CMB and faint galaxies 
as well as CMB-galaxy correlations. 



APPENDIX A: SYNCHRONOUS GAUGE 

Most public Boltzmann codes are written in the syn- 
chronous gauge (e.g. [lB,[i3])- For completeness, we out- 
line here an implementation of PPF that is appropriate 
for such codes. 

The synchronous gauge is defined by the conditions 

A = B = 0, 
1 



Vt 

hL = &Hl 



^ ' 



T = a- 



drjaA + cia 



L = - j dr]{B + kT) + C2 



(Al) 



which leave two constants Ci and C2 to be specified by 
the initial choice of time slicing and threading (e.g. the 
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rest frame of the dark matter). On the other hand the 
gauge transformation from a specific synchronous gauge 
choice to the comoving gauge used in the main paper is 
completely defined 



T 
L 



6 



hL), 



(A2) 



where the superscript s denotes synchronous gauge quan- 
tities. Note further that the comoving gauge variables 
depend only on L' and not L since vanishes in the 
comoving gauge. 

The gauge transformation equations in the main text 
allow for the evolution of the dark energy parameter F to 
be expressed in terms of the synchronous gauge variables. 
We begin by expressing the comoving gauge density fluc- 
tuations as 



^iPi = 5tP^ - p[{vT/kH) ■ 



(A3) 



This allows us to obtain $_ in terms of (5|, w| and 
r. The sources for T must also be expressed in terms of 
synchronous gauge quantities. In addition to $_ and the 
gauge invariant tt^ there is Vt- Starting with the gauge 
relation (123]) 



Vt ^ kH{<^ - : 



(A4) 



we can then use the synchronous to comoving gauge re- 
lation 



C = -m 



Iff 



and the metric relation implied by Eqn. (|32p 



(A5) 



(A6) 



H 



to obtain 



^ = ^ + + !)<!>- -^PTn. + ,., (A7) 



47rG 



Note that 

Vl = - + 4 , (A9) 

and is given by Eqn. (|40|) in terms of F, and Vt- 
Since Pe^e is gauge invariant, the remaining component 
is the pressure fluctuation obtained from 



6pl = ^Pe+P'e 



(AlO) 



The third and fourth Einstein equations fT^]) give 
1. 



— in' 
AttG 



6 

{PT 



h'L) 



PT)-r- 

kn 



(Pe 



V 



(All) 



and 



hi + 



SttG 



[StPt + iSpfp + 51pt + iSpl 



(A12) 



Since the pressure perturbation (5pg is cumbersome to 
evaluate, one can alternately replace equation (jA12p with 
the first Einstein equation ((T3)) 



h'L = 2cKk%rjT 



SttG 



(S^p^T + S^ePe): 



(A13) 



which involves w| and SePl only. The latter swtem is the 
set of Einstein equations solved by C AMB [l3| . 

The only drawback of using the latter system is that 
the matter density and metric depend to leading or- 
der on the velocities of the effective (or real) dark energy 
through TjT even at kn S> 1 since the synchronous ini- 
tial conditions typically set the dark matter velocity to 
zero. For example in ordinary smooth dark energy mod- 
els, dark energy perturbations cannot be consistently set 
to zero in this set of equations when solving for the mat- 
ter power spectrum whereas they may be in the former. 



which now gives the evolution equation (|37p for F in 
terms of synchronous variables only. 

The conservation equations for the separate compo- 
nents of real matter and radiation are of course unaltered 
but require the evolution of /i^ and rjT from the Einstein 
equations. These are modified to include the stress en- 
ergy of the effective dark energy. We can use the fluid 
correspondence and gauge relations to obtain 
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